A polytropic model of quark stars 



X. Y. Lai and R. X. Xu^ 



School of Physics and State Key Laboratory of Nuclear Physics and Technology, 
Peking University, Beijing 100871, China 



Abstract 

A polytropic quark star model is suggested in order to establish a general frame- 
work in which theoretical quark star models could be tested by observations. The 
key difference between polytropic quark stars and the polytropic model studied pre- 
viously for normal (i.e., non-quarkian) stars is related to two issues: (i) a constant 
term representing the contribution of vacuum energy may be added in the energy 
density and the pressure for a quark star, but not for a normal star; (ii) the quark 
star models with non-vanishing density at the stellar surface are not avoidable due 
to the strong interaction between quarks. The first one implies that the vacuum 
inside a quark star is different from that outside, while the second one is relevant 
to the effect of color confinement. The polytropic equations of state are stiffer than 
that derived in conventional realistic models (e.g., the bag model) for quark mat- 
ter, and pulsar-like stars calculated with a polytropic equation of state could then 
have high maximum masses (> 2Mq). Quark stars can also be very low massive, 
and be still gravitationally stable even if the polytropic index, n, is greater than 3. 
All these would result in different mass-radius relations, which could be tested by 
observations. In addition, substantial strain energy would develop in a solid quark 
star during its accretion/spindown phase, and could be high enough to take a star- 
quake. The energy released during star-quakes could be as high as ~ 10^'' ergs if 
the tangential pressure is ~ 10^^ higher than the radial one. 
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1 Introduction 



It depends on the state of matter at supra-nuclear density to model pulsar's 
structure, which is unfortunately not certain due to the difficulties in physics 
although some efforts have been made for understanding the behavior of quan- 
tum chromo-dynamics (QCD) at high density. Of particular interest is whether 
the density in such compact stars could be high enough to result in unconfined 
quarks (quark matter). Stars composed of quarks (and possible gluons) as the 
dominant degrees of freedom are called quark stars, and there is possible ob- 
servation^ evid^cethat_^^ stars could be quark stars (see reviews, 
e.g., IWebei] . l2005l : Ix^ . l2008al Jbll But it is still a problem to model a realistic 
quark star for our lack of knowledge about the real state of quark matter. 



The study of cold quark matter opens a unique window to connect three 
active fields: particle physics, condensed niatter physics, and astrophysics. 
Many possible states (see, e.g., lAlford et al.l . l2008l ) of cold quark matter are 



proposed in effective QCD models as well as in phenomenological models 
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uark matter could be in a sohd state ( Xu 



Mannarelli. Rajagopal &: Sharmal . 120071 ) 



Solid relativistic stars are challenging astrophysicists since the stelar matter 
can not be well approximately by a perfect fiuid and the conventional TOV 
(Tolman-Oppenheimer-Volkov) equation is thus not applicable. Nevertheless, 
in case of static and spherically symmetric gravity, the equilibrium equation 
could be similar to the TOV equation, by introducing a deviation between 



radial and tangential pressures (see, e.g., IXu et al.l . l2006l ). However, one has 
also to know the radial pressure, P, as a function of density p (and possible 
other parameters) in order to model a quark star in a solid state. 



No realistic relation of P{p) is available since no cold quark matter has been 
discovered experimentally and/or observationally with certainty, although many 
modeled relations between P and p are proposed in the literatures. Among the 
relations, a class of linear equations of state, P = nip—p'Y is currently focused , 
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with two free parameters k and p' (see, e.g.. lZdunikl . l2000l : lSharma fc Maharaj 
20071 ). in the framework of the bag model. Both relations derived in the bag 



model and in the density-dependent quark model can be regarded as special 
cases of the linear relation of P{p)- Whatsoever, the linear equation c ould 



not b e adequate if possible quark-clustering occurs in cold quark matter (IXu 



20031 ). Such matter with clustered quarks could be in a fiuid state at high tem- 
perature but in a solid state at sufficient low temperature. It should be worth 
noting that the interaction between quarks in a fireball with quarks and glu- 
ons is still very strong (i.e., the strongly coupled quark-gluon plasma. IShuryakl . 
20061 ). according to recent achievements of relativistic heavy ion collision ex- 
periments. Such a strong coupling may naturally render quarks grouped in 
clusters. How can then one state a reliable P-p relation in order to establish a 
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framework in which theoretical stellar models could be tested by observations 
if the quark-clustering effect is included? 



Astronomers faced a similar problem when trying to model non-quarkian nor- 
mal stars (e.g., main sequent stars and white dwarfs). A polytropic model 
with equation of state, P = Kp^ {V = 1 -|- had been extensively stud- 



ied previously (see, e.g.. IChandrasekharl . 119391 ) for main sequent stars as well 



as white dwarfs under the Newt onian gravity. This model has also been ex - 



tended under general relativity (iTooperl . Il965l : IShapiro fc Teukolskyl . Il983l ). 
The polytropic models are valuable because they could help us to model stars 
composed of realistic matter, such as ideal gas, photon gas, and degenerate 
fermi gas. Would it be possible for us to do a parallel investigation for quark 
stars? As discussed previously, QCD, which is still developing in low-energy 
regime, should be involved to describe cold quark matter and the equation 
of state of quark stars is then uncertain up-to-now. Furthermore, it might be 
problematic to calculate th e state of rnacrosc opic quark stars in QCD because 
QCD is till a local theory (iQiu &: Xul . l2006l ). Nonetheless, if quarks are clus- 
tered in quark stars where quarks are coupled strongly, the state of cold quark 
matter might be approximated phenomenologically by polytropic equations of 
state, since one may draw naively an analogy between the clusters in quark 
matter and the nuclei in normal matter. 



Some authors had applied the polytropic equations of sta te to model hy- 



brid stars wh i ch ha. ve cores composed of unconfined quarks (ILin et al.l . 12006 



Zdunik et al.l . 120061 ). However, their polytropic equations of state are only 



used to describe normal phase (composed of normal baryons) and mixed phase 
(composed of denser baryon matter); for the quark phase inside the core, the 
equation of state is still the form of the bag model. In this paper, we will 
apply for the first time the polytropic equations of state for quark matter and 
calculate the structures of quark stars, with different polytropic indices, n. 
These stellar models could also be regarded as an extension to the quark star 
models with linear equations of state. 

We are going to model quark stars in two separated ways, (i) The vacuum 



energy inside and outside quark matter is different. As proposed by iBambi 
(120071 ). quark stars could represent an laboratory to investigate the cosmolog- 
ical constant problem if some compensate field exist to generate a constant- 
like term to compensate the difference of vacuum energy, and the mass-radius 
curves of quark stars will be different from the standard case. We generalize 
Bambi's idea to discuss generic polytropic equations of state, (ii) The vacuum 
inside and outside of a quark star is assumed the same, i.e., quark stars have 
no QCD vacuum energy. In both key difference between polytropic 

quark star and normal star models lies on the surface density psur (Psur > 
for the former but psur = for the latter), since a quark star could be bound 
not only by gravity but also by additional strong interaction due to the strong 
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confinement between quarks. Analogously, a quark stars without QCD vac- 
uum energy could be similar to an asteroid with a sharp surface where the 
density is also none-zero. The non-zero surface density is still natural in the 
case with the linear equation of state, where the binding effect is represented 
by the bag constant, B (and then psnr = 45). 



The stability of a polytropic star depends also on the surface density. It is 
well known that a normal star with zero surface density should be unstable if 
n > 3 in the Newtonian gravity, and t hat the case of n = 3 is still unstable in 
general relativity (see details in, e.g., IShapiro fc Teukolskyl . Il983l ). However, 
in the models we will demonstrate, a quark star could still be stable even if 
n > 3. 



The structures o f related compact stars have been studied in general relativity 
by some authors. iDurgapal fc Bannerjil (119831 ) derived an analytical expression 
of niass-radius relation for isotropic stars in general re lativity, and Herrera et 



al. (IHerrera fc Santos 1 . 119971 : iHerrera &: Barretd . |2004| ) analyzed a set of solu- 



tions to the Einstein's equations for anisotropic matter. It is very i nteresting to 
model stars w i th an i sotropic pressu re for some physical reasons (IRuderman . 



19721: ISawverl. 



1972 



Sokolovl . Il980h . Harko and Mak dHarko fc Makl . I2OO2I : 



Mak fc Harkd . l2003l ) derived an analytical expression of mass-radius relation 
for anisotropic stars in general relativity, and discussed the constraints for 
the anisotropic parameter. They had also presented an exact analytical solu- 
tion of the gravitational equations desc ribing a static spher ically symmetric 
anisotropic quark matter distribution (IMak fc Harkd . |2002| ). These authors 
did not start with an equation of state, but studied the density and pressure 
in a more general framework based on the energy-momentum tensor. It is 
worth noting that the results obtained by above authors are still parameter- 
and assumption- dependent, even for the so-called exact solutions. We will 
alternatively study the problem, with an explicit form of equation of state. 
We also consider the case of anisotropic stars and compute the gravitational 
energy re leased during quakes of solid quark stars, with a parameter range 
given by iHarko fc Makl (120021 ). Numerical results show that, if the tangen- 
tial pressure which is slightly larger than the radial one changes abruptly, the 
gravitational energy released could be high enough to power the supergiant 
flares observed from soft 7-ray repeaters. 



This paper is arranged as follows. The details of polytropic model of quark 
stars of perfect and unperfect fluids, respectively, are presented in §2 and §3. 
The numerical results are shown in §4. The paper is concluded in §5. 
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2 Stars of Perfect fluid 



2.1 Quark Stars without QCD Vacuum Energy 



If there is no difference between the vacuum inside and outside of a quark star, 
the equation of state for a quark star is the standard polytropic model, with a 
non-zero surface density, representing the strong confinement between quarks. 
In this point of view, quark stars could analogously be similar to asteroids: 
the electromagnetic force dominates over gravity in the later, while the strong 
interaction can not be negligible in the former. Consequently, both of those 
objects can have a sharp surface where the density goes down to zero in a 
negligible small scale. 



Stars of perfect fluid in general relativity were discussed by iTooperl (119651 ). 
with an equation of state, 

P = Kp/, (1) 
p = Pgc'' + nP, (2) 



where pg is the part of the mass density which satisfles a continuity equation 
and is therefore conserved throughout the motion, and F = 1 + 1/n. In the 
static case with spherically symmetry, with the space-time metric of the form, 

ds' = e'cV - r\d^^ + sm^M^^) - eMr^, (3) 



the h ydrostatic equilibrium condition is derived to be (lOppenheimer fc Volkoffl . 



19391) 



1 -2GM(r)/cV^2rfP ^ GM{r) AnG 3p_^ 



P + pc'^ 



ar 



(4) 



and 



M 



R 

J p/(? ■ Anr^dr. 



(5) 



Similar to the Lane-Emden equation of normal stars, re-scale density pg and 
radius r, as well as M(r), by 

Pg = pg,9\ r = i/A, (6) 
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(7) 



where 



AnGpgc _ Kpgc 



and pgc is the rest mass density at the center, we can then obtain 

'-^^e' + V + afe''+^ = 0, 



1 -2{n + l)av/^^^ ^eS/in+i 



1 + {n + l)ae 

v' = ^^e''{l + nae), (9) 

with the initial conditions, 

^(0) = 1, v{0) = 0. (10) 



The mass and radius of a star are evaluated at the point when the den- 
sity reaches the surface density. In general, the surface of a star is defined 
by th e position where the pr e ssure is zero, or on ly the radial pressure is 



zero (IHerrera fc Santos I . 119971 : iHerrera et al.l . |2008| ) in the anisotropic case 



where the radial pressure is different from the tangential one. In fact, the sur- 
face in our calculations may not be the real physical surface, since the pressure 
could not be zero there. However, this does not significantly affect the mass 
and radius we obtain, because a quark star has a sharp edge, and the pressure 
and density decrease to zero in a layer with thickness of few femto-meters. 



2.2 Quark Stars with QCD Vacuum Energy 
The energy-momentum tensor has the form 

'J- '^'PS'i'ticles _j_ ^vacuum (11) 

For a constant vacuum energy, the energy-momentum tensor can be written 

as 

r~ = Kg,,, (12) 

and one finds 

^vacuum /j-vacuum ^ (13) 
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vacuum /T-vacuum a A\ 

P ~ -'00 ~ \^^) 

We can infer that no matter what forms the equations of state for the particles 
are, the contribution of vacuum is of the above form. Consequently, we write 
the equation of state as 

P = Kpg"- - A, (15) 
p = PgC^ + nKpg^^^ + K. (16) 

In this case, the density at surface (where pressure is zero) should also be 
non-zero. It is worth noting that these general equations of state of Eq.(15) 
and Eq.(16) could be simplified into a few special ones: the form without 
QCD vacuum energy if A = 0, the linear (relativistic) form if pg = 0, and the 
non-relativistic form if the second term in Eq.(16) is neglected (see the next 
sub-section for more discussions). 

The hydrostatic equilibrium conditions are also determined by Eq.(4) and 
Eq.(5). But in this case, we cannot derive differential equations for density 
and radius, and can only numerically calculate the structure of a quark star 
from the center to the surface to obtain the mass and radius. 



2.3 Comparison with MIT bag model 



In the MIT bag model, quark matter are composed of massless up and down 
quarks, massive strange quarks, and electrons. Quarks are combined together 
by an extra pressure, denoted by the bag constant B, which is the vacuum 
energy density similar to the A-parameter in the pol ytropic mode l . For the 



comparison, we apply the formulae given by Alcock (lAlcock et al.l . Il986l ) to 
calculate the equation of state, with strange quark mass = 250MeV and 
the strong coupling constant as = 0.6 for indications. 

Note that the i^-parameter in §2.2 could be calculated in the way of, from 
Eq.(15), 

K = A-p;i'''^\ (17) 



where psur is the surface density, since at the surface the pressure is zero. 
This density should be determined by the behavior of the elementary strong 
interaction. In the calculation, we assume that A has the same value as the bag 
constant B of MIT bag model, and choose the surface density, psur = l-Spo, 
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Fig. 1. To compare the equations of state discussed, including the polytropic states 
with A = 80MeV (sohd hues), A = (dashed hues), the corresponding non-rela- 
tivistic case with A = 80MeV (dash-dotted lines), A = (dotted lines), and that 
derived in the MIT bag model with the mass of strange quark m<j = 250 MeV and 
the strong coupling constant Ug = 0.6 (thin lines with dots), for a given surface 
density p^m = l-5/Oo- Here and in the following figures, pQ is the nuclear saturation 
density. It is evident that polytropic equations of state are stiff. 

where po is the nuclear matter density. The i^-parameter could be larger if 
the quark self-confinement effect is included, and the value determined by 
Eq.(17) is the minimal and the maximum masses presented in Fig. 2 and 
Fig. 3 could also be higher. For the sake of simplicity, we suppose a value 
of K from Eq.(17), and assume the i^-parameter is the same for both cases 
with and without vacuum energy under a same polytropic index n. If quark- 
cluster inside a quark star are very massive, the kinetic energy density could 
be negligible compared to the rest mass energy density, i.e., the equation of 
state could be non-relativistic, and the total energy density would include only 
the rest mass energy densitjo in Eq.(16). 

The equations of state discussed above are shown in Fig.l. An obvious con- 
clusion is that the equations of state of polytropic form could be stiffer than 
that of the MIT bag model, especially in the non-relativistic case. A stiffer 



^ An example similar to this non-relativistic equation of state is of the matter in 
white dwarfs where the energy density is dominated by the rest mass of nuclei, while 
the electron gas contributes to the pressure. 
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equation of state would lead to a larger maximum mass of quark stars. 



2.4 Gravitational energy in general relativity 



The gravitational energy in general relativity was calculated by lTooped (119651 ). 
In general relativity, the integrating of space volume should be different from 
that in the Newtonian gravity due to the space-time curvature around a mas- 
sive star. The proper energy, Eq, is obtained by integrating the energy density 
over elements of proper spatial volume, 

Eq = Eog + Eofc, (18) 



where the rest energy, Eog, of the system and its microscopic kinetic energy, 
Eok, are 

R 

Eog = Mogc" = 47r I pgC^e^/^rMr, (19) 


R 

Eok = 4:11 J nPe^/VMr, (20) 




where A can be calculated by 



2GM 



e"^ = 1 - (21) 



The total energy of a star with mass M = M{R) is E = Mc , and its gravi- 
tational energy, Q, is the difference between the total energy and the proper 
energy, 

n = Mc^ - Eo. (22) 



3 Stars with an anisotropic pressure 



Fluid within inhomogeneous pressure is imperfect, and we will consider only 
the case of spherical symmetry, that the tangential and radial press ure are not 



equal. In this case the hydrostatic equilibrium condition reads (e.g.. IXu et al. 
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20061 ) 



2GM(r)/cV 



P 



r — 2erp) + 

dr 



GM{r) 4nG 



r^P = 0, 



(23) 



where e is defined hj P± = {1 + e)P, and P is the radial pressure and Pi is 
the tangential one. 

Combine the hydrostatic equihbrium condition and equation of state, one can 
calculate the structures of quark stars with and without QCD vacuum energy. 



4 Numerical results 



Based on the formulae presented in §2.1, the mass-radius relations for various 
index, n, can be calculated. We are applying the Runge-Kutta method of 
order 4 to solve the differential equations, until the density reaches the surface 
density. For the case in §2.2, we numerically calculate from the center to the 
surface and obtain then the mass and radius. 

It is worth noting that the non-zero surface density of quark stars play an 
important role in the computation. This density should be determined by 
the behavior of the elementary strong interaction, and is then an uncertain 
parameter. In the calculation as following, we choose the surface density, psur = 
1.5po; where po is the nuclear matter density. 

In calculating the gravitational energy, we numerically integrate from the cen- 
ter to the surface for both cases with and without QCD vacuum energy. 



4-1 Mass-radius relations for stars of perfect fluid 



The mass-radius curves for both cases with and without QCD vacuum energy 
are shown in Fig. 2. It is evident from the calculation that the maximum mass 
of quark star decreases as the index, n, increases. This is understandable. A 
small n means a large F, and the pressure is relatively lower for higher values 
of n. Lower pressure should certainly support a lower mass of star. 

It could have observational implications that the maximum mass of quark 
stars with polytropic equations of state are larger than that derived in con- 
ventional model. Recently the results of 19 years of Arecibo timing for two 
pulsar s in the globular c luster NGC 5904 (M5) had been reported by Freire 
et al. (iFreire et all . l2008l ). They confirmed that for one of the binary pulsars 
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Fig. 2. Mass-radius relations for different polytropic indices, n, with psnr = I-S/Oq. 
Solid lines are A = 80MeVfm~^, dashed lines are for A = 0, dash-dotted lines are for 
non-relativistic case with A = 80MeVfm~^, and dotted lines are for non-relativistic 
case with A = 0. 



(M5B) the mass is 2.08 ± O.IOM© in la, and concluded that this mass for the 
pulsar would exclude most "soft" equations of state for dense neutron matter. 
However, a quark star in the polytropic model could be more massive than 
in previously derived realistic models (e.g., the MIT bag model) because of a 
stiffer equation of state, and the maximum mass could be larger than 2Mq. It 
is worth emphasizing that the maximum mass of quark star depends on the 
value of i^-parameter. Although in this paper we derive its minimum value 
under some assumptions, its real value is still uncertain. On this point of view, 
even the pulsars of masses larger than IAMq have been observed, the case of 
n = 3 could not be ruled out because the value of K could be larger than the 
value we use. 



On the gravitational stability. A polytropic star, with a state equation of P oc 
p^, supports itself against gravity by pressure, PR^ (note: the stellar gravity 
oc M/B? oc pR). Certainly, a high pressure (and thus large F or small n) is 
necessary for a gravitationally stable star, otherwise a star could be unstable 
due to strong gravity. Actually, in the Newtonian gravity, a polytropic normal 
star (with psur = 0) is gravitationally unstable if n > 3, but the star should be 
still u nstable if n = 3 when the GR effect is included (IShapiro fc Teukolskyl . 



19831 ). 
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Fig. 3. Mass-central density relations for different polytropic indices, n, with 
Psur = l-Spo- Solid lines are A = SOMcVfm"^, dashed lines are for A = 0, dash-dot- 
ted lines are for non-relativistic case with A = 80MeVfm~^, and dotted lines are for 
non-relativistic case with A = 0. Stars with central densities greater than that of 
stars with maximum masses are gravitationally unstable. 

A polytropic quark star with non-zero surface density or with QCD vacuum 
energy, however, can still be gravitationally stable even if n > 3. A quark 
star with much low mass could be self-bound dominantly, and the gravity 
is negligible (thus not being gravitationally unstable). As the stellar mass 
increases, the gravitational effect becomes more and more significant, and 
finally the star could be gravitationally unstable when the mass increases 
beyond the maximum mass. In order to see the central density-dependence of 
stability, the calculated mass-central density curves are shown in Fig. 3. 

Sound speed. From the P — p relation in §2.1, we can derive the ratio of sound 
speed to speed of light, which cannot be greater than 1, 



c 



dP 



n + 1 P 



n 



p + P 



< 1. 



(24) 



Similarly, also from the relation given in §2.2 we can come to 

C-^r < 1. (25) 

^ dp n p + P - ^ ^ 
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Both of the two inequations lead to 



'1 - r?)Kp'g' < nc^ 



(26) 



The equation above holds if > 1, that means that the causality keeps for 
n>l. 



4-2 Gravitational energy released during a star quake 



Based on various manife stations of pulsar-like stars, a solid state of cold quark 
matter was conjectured (IXul . l2003l ). A solid stellar object would inevitably re- 
sult in starquakes when strain energy develops to a critical value, and a huge of 
gravitational and elastic energies would then be released. One way to accumu- 
late both shear and bulk forces in a solid quark star is during an accretion pro- 
cess: strain develops remarkably in massive stars, for which the gravitational 
effect is not negligible. A solid star could additionally support the accreted 
matter against gravity by these forces, unless the forces become so strong 
that a star-quake occurs. This is the so-c alled AIQ (Accr e tion - Induc ed star- 
Quake) mechanism proposed previously (IXu et al.l . l2006l : IXul . 120071 ) . which 
might be responsible to the bursts (even the supergiant flares) and glitches 
ob served in soft 7-ra y repeaters/anomalous X-ray pulsars (see a recent review 
by iMereghettil . l2008h . 



How to calculate the energy released during an AIQ? Theoretically, anisot ropic 



fluid stars could be introduced for, e.g., the presence of typ e 3A superfiuid (iKippenhahn fc Weigert 
19901), different kin d of phase transitions (ISokolovl . Il980 i). and pion condensa- 



tion (jSawyerl . ll972l ). For the general relativistic configurations, when the inter- 
actions betwee n particles could b e treated relativistically, the fluid could also 
be anisotropic (jRudermaru . 119721 ). Previous theoretical results for anisotropic 
fluid, in a simple case with spherical symmetry, could still be adaptable to 
estimate the AIQ-released energy of a solid star. 



Base d on the analyti c al stu dy of anisotropic matter, Harko and Mak (IHarko fc Mak 



20021 : iMak &: Harkd . l2003l ) discussed the constraints for the anisotropic pa- 
rameter defined in §3, e, and presented an exact analytical solution for the 
gravitational equations of a stati c spherically symmetric anisotropic quark 
matter star (IMak fc Harkd . |2002| ). It is found that the e- value could be as 
high as 10~^. Applying the formulae in §2.4 and §3, we can obtain the grav- 
itational energy released during a star-quakes in both case with and without 
QCD vacuum energy, for example, with n = 1. The gravitational energy re- 
leased for quark stars of line ar equation of st ate with two different vacuum 
energy has been calculated in IXu et al.l (120061 ). We approximate that the rest 
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Fig. 4. The gravitational energy difference between stars with and without 
anisotropic pressures, which may be released during sequential star quakes. Solid 
lines are A = 80MeVfm~^, dashed lines are for A = 0, dash-dotted lines are for 
non-relativistic case with A = 80MeVfm~^, and dotted lines are for non-relativistic 
case with A = 0. 

energy (corresponds to the total baryon mass) Egg does not change during a 
star quake, Eog — Mqc^, since the released energy is much smaller than E^g. 

The gravitational energy difference between stars with e ^ and with e = 
are shown in Fig. 4. Three supergiant flares from soft 7-ray repeaters have 
been observed, with released photon energy being order of ~ 10^^ ergs. Our 
numerical results imply that for all the parameters we chosen, the released 
energy could be as high as the observed. 



5 Conclusions and Discussions 



Because of the difficulty to obtain a realistic state equation of cold quark 
matter at a few nuclear densities, we have tried to apply polytropic equations 
of state to model quark stars in this paper. The polytropic equations of state to 
quark stars are studied in two separated cases: the vacuum inside and outside 
quark matter is the same or not. In addition, the quark-clustering could lead 
to the non-relativistic equation of state. The differences between the those 
cases could be significant in the mass-radius relations, and may be tested by 
observations. It could additionally provide a way to probe the properties of 
QCD vacuum. 

The polytropic equations of state are stiffer than that derived in previous 
realistic models, so they could lead to more massive quark stars with masses 
> 2Mq. Consequently, even when some massive pulsars have been observed, 
it still can not rule out the possibility that pulsar-like stars are quark stars. 
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Though a normal star with zero surface density can only be gravitationally 
stable if the polytropic index n < 3, a quark star with non-zero surface density 
could still be stable even if n > 3. A solid quark star may break if its strain 
energy develops to a critical value, and we calculate the gravitational energy 
released during quakes and find that the energy could be as high as ~ 10^'' 
ergs if the anisotropic parameter, e, could be order of 10~^. Such a huge 
of energy would be liberated during an AIQ (accretion-induced star-quake) 
process, to be probably responsible to the bursts and glitches observed in soft 
7-ray repeaters/anomalous X-ray pulsars. The general relativity effect has 
been included to simulate the polytropic quark stars. 



The nature of pulsars is unfortunately still a matter of controversy, even more 
than 40 years after the discovery of pulsar. Although quark stars cannot be 
ruled out, both theoretically and observationally, and pulsars are potential 
idea laboratories to study the elementary strong interaction, we are lacking 
a general framework in which theoretical stellar models could be tested by 
observations. Polytropic quark star model is the one we try to establish. Future 
advanced observations may help to constrain the uncertain parameters, e.g., 
the polytropic index n, the coefficient K, the surface density psur, and even 
the vacuum energy A. 



One of the daunting challenges nowadays is to understand the fundamental 
strong interaction between quarks, especially the QCD in the low-energy limit, 
since the coupling is asymptotically free in the limit of high-energy. The state 
of cold matter at a few nuclear densities is still an unsolved problem in the low- 
energy QCD. In effective QCD models, BCS-type quark pairing was proposed 
to form at a Fermi surface of cold quark matter, and the shear moduli of the 
rigid crystalline color super- conducting qua rk matter could be 20 to 1000 times 
large r than those of neutron star crusts (iMannarelli. Rajagopal fc Sharmal . 



20071 ). However, quark cluste rs are phenomenologically suggested to form in 



cold quark matter (IXul . l2003l ). The state of such cold quark matter might be 
approximated by polytropic equations of state since one may draw naively 
an analogy between the clusters in quark matter and the nucleus in normal 
matter. Certainly, it would be very interesting to observationally distinguish 
between those two kinds of solid quark matter. 
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